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ABSTRACT 

Using a new numerical model for cosmic chemical evolution, we study the influence 
of hypernova feedback on the star formation and metal enrichment history of the uni- 
verse. For assumptions which produce plausible results in idealized collapse models of 
individual galaxies, our cosmological simulations of the standard ACDM cosmology 
show a peak of the cosmic star formation rate at 2: ~ 4, with ^ 10% of the baryons 
turning into stars. We find that the majority of stars in present-day massive galaxies 
formed in much smaller galaxies at high redshifts, giving them a mean stellar age as 
old as 10 Gyr, despite their late assembly times. The hypernova feedback drives galac- 
tic outflows efficiently in low mass galaxies, and these winds eject heavy elements into 
the intergalactic medium. The ejected baryon fraction is larger for less massive galax- 
ies, correlates well with stellar metallicity, and amounts to ~ 20% of all baryons in 
total. The resulting enrichment history is broadly consistent with the observed abun- 
dances of Lyman break galaxies, of damped Lyman a systems, and of the intergalactic 
medium. The metallicity of the cold gas in galaxies increases with galaxy mass, which 
is comparable to observations with a significant scatter. The stellar mass-metallicity 
relation of the observed galaxy population is well reproduced by the simulation model 
as a result of the mass-dependent galactic winds. However, star formation does not 
terminate in massive galaxies at late times in our model, and too few dwarf galaxies 
are still forming stars. These problems may be due to a lack of resolution, to inap- 
propriate modelling of supernova feedback, or to a neglect of other feedback processes 
such as active galactic nuclei. 
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1 INTRODUCTION 

While the evolution of the dark matter in the standard 
ACDM cosmology is reasonably well understood, the evo- 
lution of the baryonic component is much less certain be- 
cause of the complexity of the relevant physical processes, 
such as star formation and feedback. One approach to study 
the dynamics of baryons i n galaxy formation is to use semi- 
analytic techniques (e .g., iKauffmann. White fc Guiderdonil 
ll993l : ICole et al.lll994 ) that combine the growth of dark mat- 
ter halos with simplified parametorizations of the baryonic 
physics. While this allows the construction of theoretical 
models for galaxy formation and chemical enrichment with 
a large dynamic range, the validity of the assumptions must 
be ultimately justified by observations and by more detailed 
theoretical work. We here adopt the other common approach 
to study galaxy formation, which is to use direct hydrody- 
namical simulations. Complementary to the semi-analytic 
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models, they treat the dynamics of baryons in a much more 
detailed fashion. 

Many simulation codes have been developed for 
studying galaxy formati on and evo l ution, not only of 
isolated systems (e.g., Kat^ 199?, Mihos & Hornquis 
J,989, Stcinmctz & Miillcr 1994, Nakasato & Nomoto 200._ 
Ikawata fc GibsonI |2003| . iKobavashil |200j) but also for 
cosmological s imulations of individual galaxies (e.g., 
iNavarro fc Whi te 1994) or of the galajcy population a s a 
whole l|Cen fc Ostrikeni 19991 . [Springel fc Hernquistll2003l ). A 
robust result has been that with the commonly employed, 
schematic star formation criteria alone, i.e., converging gas 
fiow, local Jeans instability, and short cooling times, the pre- 
dicted star formation rates (SFRs) are higher than what is 
compatible with the observed luminosity density. Thus feed- 
back mechanisms are in general invoked to reheat gas and 
suppress star formation. We include both supernova and 
hypernova feedback in our hydrodynamical model in this 
paper. Supernovae inject both thermal energy and heavy 
elements into the interstellar medium, and by means of 
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supernovae-driven galactic winds, some of these metals can 
escape from galaxies into the intergalactic medium. 

There is observational evidence for galactic winds 
star fo rming galaxies both lo cally and a t high 



redshifts llHeckman et al.1 I2q00l; iPettini et all [20011; 



Martin. Kobulnickv fc HeckmanI l2002l : iMazzotta et al.l 
20021 : lOhvama et al.ll2002r ). The existence of heavy elements 



in the intracluster medium (ICM) and intergalactic medium 
(IGM) also requires such an ejection mechanism. The 
mass-metallicity relation and also the color-magnitude 
relation of early-type galaxies have been explained as 
the t he result of loss of heavy elements in galactic winds 
(e.g., iLarsonl [l973 . lArimoto fc Yoshiil [l987l '). The relation 
between the iron m ass in the ICM and the luminosity of 
early-type galaxies IjArnaud et al.l 119921 ) may also suggest 
metal ejection from early-type galaxies. 

Since the presenc e of heavy elements can inc rease gas 
cooling substantially (jSutherland fc Dopital 119931 ). consis- 
tent simulations of galaxy formation must include a treat- 
ment of chemical enrichment as well as energy feedback. 
However, especially with respect to chemical enrichment, 
most existing hydrodynamical models are too simplistic to 
be compared in detail with the real Universe. There are only 
a few studies that simulate the detailed chemical enrich- 
me nt by both Type II and la supernovae (SNe II and SNe 



la) l|Raiteri. Villata fc Navarro! 



199a iMosconi et al.l 



200 ll. 



1996, Carr aro. Lia fc Chiosil 
Nakasato fc Nomotol 12001 



Kawata fc GibsonI 12003 , iTornatore et al.l |2004 iKobavashil 



200J). Different types of supernovae produce different heavy 



elements on different timescales. Although the most abun- 
dant heavy element is oxygen which accounts for half of the 
solar metallicity, this element is hardly detected in observa- 
tions of stars. Iron, which is another fundamental element 
observed in the damped Lyman a (DLA) systems and the 
ICM, is mainly produced by SNe la. Recently, it has been 
confirmed that supernovae with ten times the standard ki- 
netic energy exist and produce a certain amount of iron, 
not only by the ob servation of individual supernovae (e.g., 
iNomoto e t al. 2002) but also by the ele mental abundance ra- 
tios of Milky Way stars (JKobavashi et al. 2006). Even with 
the Salpeter initial mass function (IMF), three times more 
energy is ejected on average when such hypernovae (HNe) 
are included, which therefore increase the overall strength 
of supernova feedback substantially. 

In the context of hierarchal galaxy formation, a number 
of important questions should be explained by simulation 
models, including (1) the exist ence of massive ga laxies at 
high redshifts (e.g., Forstcr Sc hreiber et al.ll2004 ). (2) the 
old stellar populations of massive elliptical galaxies (e.g., 
iBower. Lucev fc Ellislll992l). (3) the "down-sizing" effect o f 



star formation (e.g.. lCowie et al.lll996l . lKodama et al.ll2004r ). 
and (4) the number of spiral galaxies and the detailed mor- 
phological mix of galaxies. With our simulations, we can 
only address questions (2) and (3) because of limitations in 
volume and resolution. Note that while more massive halos 
tend to form at later times in hierarchal clustering, this does 
not necessarily imply that massive galaxies have younger 
stars than low-mass galaxies. This is because the present-day 
massive galaxies may have formed by the merging of smaller 
galaxies, i.e. the assembly time of these systems needs not 
to be the same as the formation time of their stars. To study 
this question, we will directly analyze when and where stars 



form in our simulations and we discuss the ages of stellar 
populations of massive galaxies formed in our models. 

It is not clear that the IMF is universal, and simulations 
of the formation of Population III stars appear to suggest a 
flatter IMF or a l ack of low-mass stars in low-metallicity en- 
vironments (e.g., lBromm fc Larsonll2004 ). This uncertainty 
means that the SFR and the stellar mass cannot be uniquely 
obtained from the observed luminosity. However, since a flat- 
ter IMF produces more heavy elements, the metal enrich- 
ment history of galaxies can put constraints on their star 
formation history and IMF. In fact, the observed weak chem- 
ical evolution of the IGM (Songaila 2001, Schayc ct al. 200j) 
and of DLA systems (jPettini et al.l 1 19971 . iProchaska et al.l 
[2003) rnay conflict with a flatter IMF. A detailed treatment 
of chemical enrichment is required in cosmological simula- 
tions. To this end we will analyze where and how heavy 
elements are produced and distributed in the gas phase and 
in stars, using a Salpeter IMF for definiteness. 

In this paper, we simulate the evolution of gas and stel- 
lar systems from cosmological initial conditions, following 
chemical enrichment fro m SNe II, SNe la, and HNe. We use 
the GADGET-2 code by Springel (200 5), in which we intro- 
duce metal-dependent cooling rates (ISutherland fc Dopital 
I19931 ) and the chemical enrichment model of Kobayashi 
(2004, hereafter K04). The feedback scheme is slightly mod- 
ified from K04, as described in Section 2. In Section 3, as a 
test of our code, we simulate galactic winds from dwarf spi- 
ral galaxies and examine their parameter dependencies. Our 
cosmological simulations are described in Section 4, where 
we also discuss their predictions for the star formation his- 
tory, galactic winds, and the chemical enrichment. Section 5 
gives our conclusions. 



2 SIMULATION MODEL 

To construct a self-consistent three-dimensional chemody- 
namical model, we have introduced various physical pro- 
cesses associated with the formation and evolution of stel- 
lar s^^stems into the parallel tree-SPH code GA DGET- 
2 (|Springell [2005I . ISpringel. Yoshida fc White! l2000l ). which 
is fully adaptive with individual smoothing lengths and 
timest eps and uses an entropy-co nserving formulation of 
SPH (ISpringel fc Hernguistl I2OO2! ). We include radiative 
cooling, and photo-he ating by a uniform and evolving UV 
background radiation (JHaardt fc Madaulll996l ). Star forma- 
tion, feedback from Type II and la supernovae (SNe II and 
SNe la), stellar winds (SWs), and chemical enrichment are 
modelled as well, with an implementation that is close to 
the one in K04. In brief, the characteristics of this part of 
our code may be summarized as follows. 

(i) Radiative cooling is computed using a metallicity- 
dependent cooling function. For primordial gas ([Fe/H] 
< —5), we compute the cooling rates using the two- 
bo dj;_grocesses_of_jI_and_jIe, and f ree-free emission, as 
in iKatz. Weinberg fc Hernguistl l|l996! ) . For metal enriched 
gas ([Fe/H] ^ —5), we use a metallicity-dependent cool- 
ing function computed with the MAPPINGS HI software 
ijSutherland fc Dopital ll993l V In this cooling function, the 
elemental abundance ratios are set to be constant for given 
[Fe/H] according to the relations found in the solar neigh- 
bourhood. [0/Fe]= 0.5 for Galactic halo stars for [Fc/H] 
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Figure 1. The time evolution of the formation of an isolated disk galaxy in a halo of mass of 10 ft^^M©. The black points show star 
particles, while the gas particles are colour-coded according to their temperature. Each panel is 20 kpc on a side. The upper row shows 
face-on projections, the lower row gives edge-on views. 



^ —1, and solar values for [Fe/H] ^ 0. We interpolate be- 
tween these values for —1 < [Fe/H] < 0. 

(ii) Our star formation criteria are the same as in lKata 
l|l992h : (1) converging flow, (V • v)i < 0; (2) rapid cool- 
ing, tcooi < tdyn, and (3) Jeans unstable gas, idyn < isound- 
The star formation timescale is taken to be proportional 



to the dynamical timescale (fsf 



-*dy 



where 



star formation t imescale parameter which we set to 0.1 
l|Kobavashil |2005| ) . For every timestep At, we draw a ran- 
dom number p between and 1, and provided it fulfils 

a fractional part of the mass 



P 



.{ 



exp 



At 



of the gas particle turns into a new star particle. The mass 
of the new star particle is given as m*,o — m^fi/N, with 
A^, = 2. Note that an individual star particle has a typical 
mass of ~ lO^M©, i.e. it does not represent a single star but 
an association of many stars. The masses of the stars asso- 
ciated with each star particle are distributed according to 
an initial mass function (IMF). We adopt a power-law IMF, 
(?i(m) oc m~^ (the slope x = 1.35 gives the Salpeter IMF), 
which is assumed to be independent of time and metallicity. 
We limit the IMF to the mass range O.IMq < m sJ 12OAf0. 

(iii) We do not adopt the instantaneous recycling approx- 
imation when accounting for feedback. Instead, we treat 
star particles as evolving stellar populations which eject 
thermal energy Ee, gas mass Em, and heavy elements E^. 
from SWs, SNe II, and SNe la, as a function of time. The 
release of energy and heavy elements is distributed to a con- 
stant number A'^fb of surrounding gas particles. We will show 
results for two cases: In our 'SN feedback' model, the SN II 
energy and yields are the same as adopted in K04. In the 
'HN feedback' model, we assume instead a constant hyper- 
nova fraction ehn = 0.5 for progenitor masses M ^ 20 Mq, 
and adopt the mass-en ergy relation; 30, 20, 10 , 10 x 10^^ erg 
for 40, 30, 25, 20 Mq (iKoba vashi et al.ll2006t l. The ejected 
energy of each SW, SN II, and SN la are ~ 0.2 x lO'^^ erg 
depending on metallicity, 1 x 10^^ erg, and 1.3 x 10^^ erg, 
respectively. 

We distribute this feedback energy in purely thermal form, 
although a fraction of it could, in principle, be distributed 



in kineti c form as a velocity per turbation to the gas parti- 
cles (see lNavarro fc Whitelll993h . The energy increase of a 
gas particle i due to a star particle j is calculated in each 
timestep Atj as 



'df 



/ dtE,At,Z,)Wiu,)/J2wirj 



(1) 



The feedback neighbour search needs to be done twice in 
order to ensure proper mass and energy conservation; first 
to compute the sum of weights for the normalization, and 
a second time for the actual distribution. Roughly speak- 
ing, since we assume constant A^fb, the energy increase 
per mass is almost constant and independent of resolution; 
Ee/mg/NFB -- m,,o7lsNl0^^erg/mg/AfFB oc 1/7V./7V"fb. 
(Naively, the mass of gas particles are decreasing because of 
star formation.) The feedback radius that was assumed to 
be constant in the GRAPE-SPH code of K04 is now variable 
depending on the local gas density. With higher resolution, 
constant A^fb gives smaller feedback radius. 

For the metals, the mass- and metallicity-dependent nu- 
cleosynthesis yields of SNe II and SNe la are taken from 
Kobayashi et al. (2006) and Nomoto et al. (1997b), respec- 
tively. The progenitor mass ranges of SWs and SNe II are 
8 - 120 Mq and 8 - 50Mq, respectively. For SNe la, we 
adop t the single deg enerate scenario with the metallicity ef- 
fect (jKobavashi et al . 1998, 2000), where the progenitors are 
Chandrasekhar white dwarfs (WDs) with an initial mass of 
3 — 8 Mq, and the lifetimes are determined from the lifetimes 
of the secondary stars with 0.9 — 1.5 Mq and 1.8 — 2.6 Mq 
for the red-giant (RG) and main-sequence (MS) systems, 
respectively. The fractions of WDs that eventually produce 
SNe la are adjusted to match the chemical evolution con- 
straints of the Milky Way Galaxy as [foRc = 0.02, bus ~ 
0.04]. 

(iv) The photometric evolution of a star particle is 
identical to the evolution of a simple stellar p o pulat ion. 
Spectra f\ are taken from iKodama fc Arimotd (|l997l ) as 
a function of age t and metallicity Z. 
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Table 1. Properties of different simulations of isolated disk galaxy formation, at time 4 = 13 Gyr. In all runs, we adopted c = 0.1. The 
total mass and gravitational softening length are given in units [/i^^IO^^Mq] and [h~^ kpc], respectively. The 'wind gas' is here defined 
as the gas particles that are outside r > 2r2oo- /w and /z,w denote the wind fraction relative to the total baryon mass and the ejected 
metal fraction relative to the total amount of metals produced, respectively. 



3 ISOLATED DISKS 

As a test of our code, we show the evolution of an iso- 
lated disk galaxy in a static dark matter potential. W e con- 
struct an NFW-halo (|Navarro. Frenk fc Whit^ligQeD with 
10% gas content in virial equilibrium for total masses of 
Mvir = lO^°/i"^M0 and 10^^ h'^MQ. The initial angular 
momentum corresponds to a spin parameter A = 0.1 and is 
distributed within the halo under the assumption that the 
specific angular mom enta of spherical shells are all aligned 
l|Bullock etal.l 120011 ). For the IO^^/i'^Mq model, we work 
with a numerical resolution of A^gas ~ 10000, 40000, and 
160000 particles, corresponding to a gas particle mass of 
1 X 10^ 2.5 X lO"*, and 6.25 x 10^ /i"^MQ, respectively. The 
adopted gravitational softening length is given in Table [l] 

Figure [T] shows the time evolution of gas and star par- 
ticles in the 10^" /i~^Mq halo for the case with SN feed- 
back (A'' — 160000 resolution). After the start of the sim- 
ulation, the gas is allowed to cool radiatively, making the 
gas in the centre of the halo quickly lose its pressure sup- 
port such that it settles into a rotationally supported disk 
that grows from inside out. The rapid initial growth of the 
gas disk induces an initial starburst. Figure [5] compares the 
star formation rates (SFRs) for the models without feed- 
back (dotted line), with feedback (dashed line), and with 
metal-dependent cooling (solid line). The SFR is peaked 
at t ~ 0.07 Gyr with maxima of 0.25, 0.45, and 0.8 M© 
yr"'^, respectively. Since the shortest lifetime of SNe II is 
~ 4 X 10® yr, star formation is affected by the feedback dur- 
ing this phase. Although the adopted star formation scheme 
is different, our model without feedback is similar to that of 



ISpringel fc Hernguistl (120031 ) . With metal-dependent cool- 
ing, star formation is more efficient with a twice larger peak 
SFR. 

As soon as stars form in the centre, thermal energy 
of SNe II is injected into the inter-stellar medium. The 
energy from a star particle is isotropically distributed to 
the surrounding gas particles in our model. Because of 
the dense disk formed and the non-isotropic infall of gas, 
the low-density hot-gas region expands in a bipolar flow 
(t '^ 0.14 Gyr) as result of the energy input. In the disk 
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Figure 2. Star formation rates in isolated halos with total mass 
10^" h~l M© (upper panel) and 10^^ h~^ Mq (lower panel). The 
dotted, dashed, and solid lines show the results with no feedback, 
feedback, and feedback plus metal-dependent cooling. 



plane, the hot gas region expands and forms a dense shell 
where stars keep on forming. The energy from these stars 
can quickly propagate to the surrounding low-density re- 
gion. After the galactic wind forms (t ~ 0.28 Gyr), the gas 
density becomes so low at the centre that star formation is 
terminated for a while. Because of radiative cooling, a part 
of the ejected gas however returns, again settling in the disk 
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Figure 3. The time evolution of the wind mass in 10^ /i~ Mq 
disks with different resolution and feedback parameters for 
the model with feedback and metal-dependent cooling; A^ = 
40000, iVpB = 72 (soUd line), A^fb = 36 (dashed Une), 144 (dotted 
line), N = 10000 (dot-dashed line), and N = 160000 (dash-dot- 
dot-dot line). In the lower panel, the wind mass is normalized 
by the stellar mass at each time, which corresponds to the wind 
efRciency. 



where it fuels new star formation. This secondary star for- 
mation is not as strong as the initial starburst, and takes 
place in a more continuous fashion. Although some small 
bubbles are forming in the galaxy in this stage, not much 
gas is ejected by them from the disk. 

Supernovae eject not only thermal energy but also 
heavy elements into the interstellar medium. Since the cool- 
ing function depends strongly on metallicity, the star for- 
mation history is altered when metal-dependent cooling is 
accounted for. This is important especially for the secondary 
star formation when the gas has become enriched. Indeed, 
when metal-dependent cooling is included, the secondary 
star formation is much stronger than without it. Compared 
with the no feedback case, the final stellar mass decreases by 
a factor of 3 if we include the feedback, but increases back 
to the original value if we add metal-dependent cooling as 
well. In Table [ij we summarize the properties of the galaxies 
in the different cases after a time of 13 Gyr. 

For the massive halo with lO^'^ /i~^Mq (the lower panel 
of Fig. [2]), a galactic wind never occurs, and star forma- 
tion takes place continuously. When we include feedback, 
the SFR shows repeated bursting behaviour. However, with 
metal-dependent cooling, star formation again takes place 
smoothly. This is because the enhancement of the star for- 
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Figure 4. Surface brightness profiles as a function of projected 
radius (left) and height (right). The lower panels show the oxygen 
abundance gradients. 



mation due to the existence of heavy elements is larger than 
the suppression by thermal feedback. 

We note that these initial conditions are somewhat arti- 
ficial since galaxies are not born as pure gas disks in reality. 
We adopt these initial conditions in order to explore the ef- 
fects of feedback in a clear way, without having to deal with 
the complications of the hierarchical formation of galaxies. 
In reality, massive galaxies form through the merging of sub- 
galaxies, and galactic winds can occur from the subgalaxies 
before merging. Weather the ejected materials fall onto the 
final galaxies or not is another issue, and should depend 
on the details of the merging histories. In Section 4, we will 
adopt more realistic initial conditions appropriate for cosmo- 
logical simulations where the hierarchical formation is fully 
taken into account. 

The strength of the galactic winds depends on i) the 
feedback parameter and ii) the numerical resolution. Fig- 
ure [3] shows the time evolution of the wind mass normalized 
by the total mass (upper panel) and the stellar mass (lower 
panel) as a function of time. The wind particles are here 
defined as those gas particles that are outside r > 2r2oo- 
In our standard model, almost half of the gas is ejected by 
the wind (/w = 0.54). In our feedback scheme, the feedback 
neighbour A'^fb controls the ejected energy from one star par- 
ticle to one gas particle. Once the gas particles are heated 
to T > 3 X 10"" K where the cooling rate becomes low, they 
can avoid rapid cooling. Therefore, the feedback effect be- 
comes stronger with smaller A^fb. With A'^fb ~ 36, the wind 
mass increases to give /„ = 0.8, and the wind efficiency nor- 
malized by the stellar mass is ~ 5 times larger than for the 
A^'fb = 72 case. On the other hand, with A'fb = 144, the 
distributed energy is too low to generate a galactic wind. 
In addition, a large value of A''fb corresponds to a very 
large feedback radius (~ 10 kpc in these simulations), which 
causes unrealistically large mixing of gas. 

A comparable amount of heavy elements can be ejected 
by the wind with the feedback. However, with the metal- 
dependent cooling, the ejected metal fraction relative to the 
total amount of metals produced is as small as a few per- 
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Figure 5. The time evolution of our cosmological simulation in a periodic box 10 h^^ Mpc on a side. We sliow the projected density of 
gas (first row) , stellar V-luminosity (second row) , gas metallicity log Zg /Zq (third row) and stellar metallicity log Z, /Zq (forth row) for 
the high resolution simulation {N = 96^, A^fb = 405). 



cent for the SN feedback case. For the case with HNe, the 
ejected metal fraction is as large as /z,w = 0.25. This is 
because more metals are ejected with more energy, and en- 
riched gas can be heated enough to prevent rapid cooling. 
In this case, the wind metallicity is logZ/ZQ ~ — 1 while 
the stellar metallicity is ~ —0.4 (see Table [!}. 

In principle, our feedback scheme does not depend on 
numerical resolution because the ejected energy per unit 
mass from star particles is always the same. However, the 
star formation criteria may introduce some resolution de- 
pendence. In our disk formation simulations, the feedback 
efficiency depends on the strength of the initial star burst. 
We find that the SFR during the initial starburst is smaller 
with N = 160000 than N = 40000. This is because the Jeans 
instability criterion is not as easily satisfied in high resolu- 
tion simulations with smaller smoothing length. A^ = 10000 
is too small to simulate the galactic wind properly. As a re- 
sult, the wind mass decreases with higher resolution. How- 
ever, when we normalized the wind mass by the stellar mass, 
the present wind efficiency is the same. 

The feedback also affects the structure of the stellar 
disk. In Figure |4l we show the surface brightness profiles of 
the stars in the 10^° h~^MQ halo as functions of the pro- 



jected radius (left panels) and the height at r = 2 — 3 kpc 
(right panels). Without feedback (dotted line), the disk is 
small and thin with scale lengths of re = 0.5 and h = 0.1 kpc, 
respectively. With feedback and metal-dependent cooling 
(solid line), the disk size becomes larger (re = 1.3), while 
the thickness does not change {h = 0.13). At late times, star 
formation takes place mainly in the outer regions (r ~ 2 — 3 
kpc) in this case. Thus the surface brightness profile shows 
a bump around r ~ 2 kpc. Strong metallicity gradients are 
seen at r > 2 kpc and z > 0.2 kpc. Without metal-dependent 
cooling (long-dashed line), the disk size is smaller and the 
thickness is larger. In this case, star formation takes place 
in the centre at later times, and no metallicity gradient is 
formed. Numerical resolution does not affect these results 
(dot-dashed line). 



4 COSMOLOGICAL SIMULATIONS 

We focus on a ACDM cosmological model with parameters 
Ho = 70kms-^Mpc-\ fi™ = 0.3, Oa = 0.7, f^b = 0.04, 
n — 1, and as — 0.9. The initial conditions are set up 
in a 10/i~^Mpc cubic box with periodic boundary condi- 
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tions, employing an equal number of dark matter and gas 
particles. We work with two different resolutions; A^dm = 
A'gas = 54 and 96 . The mass of a dark matter particle 
is 4.58 X 10* Mq and 8.16 x lO'^ M© and the mass of a gas 
particle is 7.05 x lO'^M© and 1.25 x lO'^M© for these two 
resolution, respectively. The adopted gravitational soften- 
ing lengths are given in Table [2] For gas particles we adopt 
half the value adopted for dark matter. With several low- 
resolution simulations we examine dependencies on our two 
most important model parameters, the feedback neighbour 
number A'fb and the star formation timescale c*. 

In order to identify galaxies , we apply the same method 
as lSpringel fc Hernauistri|2003bl ). A friend-of-friends (FOF) 
group-finding algorithm is applied to the dark matter parti- 
cles, using a fixed comoving linking length equal to 0.2 times 
the mean inter-particle spacing of the dark matter particles. 
Then gas and star particles are associated with the near- 
est dark matter particles. We discard all groups with fewer 
than 32 dark matter particles. The intergalactic medium is 
coarsely defined as the particles that do not belong to any 
identified halo. We note that our results are insensitive to 
the detailed parameters of the FOF group finding methods. 

Because of practical limitations with respect to the 
available computational time, the box size of our simula- 
tions is quite limited. Thus very massive galaxies and galaxy 
clusters forming at low redshift are not included and the 
large-scale structure on the scale of the box-size is affected 
by non-linear evolution below z ~ 2 — 3. As a result, the 
simulated cosmic SFR in Figure |8] at z < 2 — 3 is prob- 
ably an underestimate of the cosmic mean expected for a 
larger simulation volume. However, the evolution of indi- 
vidual galaxies is expected to be less affected by the finite 
box size. In particular, the mass-metallicity relation (Figs. 
1191 and I20|) and the wind-mass relation (Fig. I16|) should be 
quite robust for low and intermediate-mass galaxies, even 
though very massive galaxies are not included. 

Figure[5]shows the time evolution of the densities of gas, 
stellar V-band luminosity, gas metallicity, and stellar metal- 
licity for the high resolution simulation (A'' = 96^,A''fb = 
405) in our 10/i~^Mpc box. Star formation takes place in a 
distributed fashion at high redshifts, and becomes most ac- 
tive around z — 3. The distribution of stars appears smooth 
at high redshifts, but concentrated at lower redshifts. Be- 
cause of the feedback, massive galaxies are surrounded by 
hot gas, and heavy elements are also distributed in the inter- 
galactic medium. Both in the gas phase and stars, metallicity 
gradients are generated in high density regions. 

Figure |S] shows the density-temperature diagram at 
z = 0, representing metallicities with colours. Metal-rich 
gas particles populate mainly the cold dense region, with a 
few particles found in the hot, low-density region. We note 
that one-tenth of particles are randomly chosen and plotted 
in this diagram, and the statistics of the low-density cold 
gas is much better. Figure [7] shows that the relation that 
the gas metallicity is higher for high density is present at all 
redshifts. This is also responsible for generating the radial 
metallicity gradients in galaxies both for gas and for stars. 



4.1 Model Comparison 

In Figure [8^, we show how feedback and metal-dependent 
cooling affect the cosmic star formation history. The cos- 




Figure 6. The density-temperature diagram at z 
coded by the metallicity of gas particles. 



0, colour- 




Figure 7. The metallicity-density relation at 2 = (solid line), 
z = 1 (dashed line), z = 2 (dotted line), 2 = 4 (dot-dashed line), 
and 2 = 5 (dot-dot-dot-dashed line). 



mic SFRs that are directly measured from the ages of 
stellar particles for the low resolution simulations (A'' = 
56'^,A^FB = 72). Supernova feedback (dot-dashed line) de- 
creases the SFR from 2; ~ 3, and the SFR is smaller by 
factor of 2 at < 2 < 2 than in the no feedback case (dotted 
line). However, metal-dependent cooling (long-dashed line) 
increases the SFR back to the no feedback case. For com- 
parison. Figure [8^ also shows observational estimates of the 
cosmic SFR density at different epochs. Note that these de- 
terminatio ns are derived from the observed luminosity den- 
sities (e.g.. lMadau et al.lll996r ). and this involves uncertain- 
ties from dust extinction and completeness, as well as from 
the IMF. Rest-frame UV observations are plotted with dust 
correction with a factor of 2.7 for 2 < 2 an d 4.7 for z > 2 
(Stcidcl ot al. 1999, Re ddv fc Steidell 120041), and dust cor- 
rected Ha measurements are taken from lSchiminovich et al.l 

liooi). 

If we include hypernova feedback (dashed line), the SFR 
starts to be suppressed from 2 ~ 6 onwards, and is overall 
smaller by a factor of 3 at < 2 < 3. The resulting SFR is 
in broad agreement with the observations that show a peak 
of log(SFR/[MQyr~'^]) ~ — 1 at 2 ~ 3. Since a hypernova 
releases ten times more energy than an ordinary supernova, 

yr) is as 
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Table 2. Mean values at present epoch z = 0. The box size and the softening length e are in [h~^ Mpc] and [h~^ kpc], respectively. 
Wind mass is calculated from tracing the orbit of gas particles, /w and /z,w denote the wind fraction relative to the total baryon mass 
and the ejected metal fraction relative to the total amount of metals produced, respectively. 



large as ~ 10^ erg, and the temperature increase of indi- 
vidual gas particles reaches ~ 10^ K. Once gas particles are 
heated to this temperature, the gas particles do not cool 
rapidly due to the comparatively low cooling rate at this 
temperature. This reduces not only the SFR but also the 
metal production rate, even though the ejected metal mass 
from a single HN is larger than from a SN. The redshift 
where the mean gas [Fe/H] reaches ~ — 2 is z ~ 3 for the 
SN feedback, but is delayed to z ~ 2 for the HN feedback. 
Since the metallicity dependence of the cooling rate clearly 
appears for [Fe/H] > — 2, metal-dependent cooling does not 
play such a prominent role in the HN feedback case, and 
enriched gas can remain hot without forming stars. 

In Figure [SJd, we show the time evolution of the stel- 
lar mass density. To avoid uncertainties in the completeness 
from the faint end, we measu re the stellar mass in the galax- 
ies with M. > 10* Mq, as in lFontana et all l|2006l ). The HN 
feedback improves the agreement wit h these observations . 
The luminosity-limited observations (jRudnick et alj |2003| . 
l200a ) are also reproduced with the luminosity-limited mea- 
surements in our simulations. 

The present stellar fraction (relative to total baryon 
density) and mean metallicities are tabulated in Table [21 
Without feedback, 25% of baryons turn into stars, which 
is too large c ompared with observational estimates (e.g., 
iFukueita fcTe cblcs 2004). With SN feedback, the stellar 
fraction reduces to 10 — 15%, which may be consistent with 
observation. Recently, the observational estimate h as been 
reduced to less than 10% (JFukueita fc Peeble^l2004l ). which 
may require larger feedback. The larger energy ejection by 
HNe could provide a solution. The present mean gas metal- 
licity is [Fe/H] ^ —1.3 for SN feedback, and is reduced to 
--^ — 1.8 for HN feedback. The mean stellar metallicity is al- 
most solar for the SN feedback, and becomes sub-solar for 
the HN feedback. 

The feedback parameter A^fb does only weakly influ- 
ence the cosmic SFRs. With A'^fb ~ 36, the SFR increases 
by a factor of 2, which is in the opposite sense from the 
isolated disk case (Section 3). With the resolution of our 
cosmological simulations, the ejected energy from one star 
particle is large, and smaller A^fb results in weaker feedback 
due to a smaller feedback radius. For the higher resolution, 
because more small dense regions are resolved, the SFRs at 



z > 4: become much larger and the SFR at 2; ~ 3 increases 
by a factor of 3. To obtain similar feedback radius, we adopt 
A'fb = 405 for the A'^ = 96^ case, which gives a similar peak 
SFR and stellar fraction to the low resolution case. 

We now discuss the dependence of our results on the 
star formation timescale parameter c*. In the one-zone 
model, the star formation timescale directly affects the star 
formation history, and it may be imagined that a smaller 
c, can easily suppress the SFR. However, in our chemody- 
namical model, the influence of c, is more subtle than that. 
With a smaller value for c*, say c, = 0.02, the peak redshift 
of the SFR is shifted towards lower redshift {z ~ 3) and the 
present SFR density is enlarged (log SFR = —1.4), but the 
total amount of stars is not reduced in any significant way. 
The reason is that if we prevent the gas from forming stars 
at some early epoch, this will only delay the transformation 
of cold gas into stars, but the amount of gas that cools and 
becomes available for star formation is not changed. If we 
pick a value as large as c, = 1, the initial star burst is so 
large that the SFR continuously decreases from z = 6. Com- 
paring with the observed cosmic SFRs, we therefore adopt 
c, — 0.1, which roughly reproduces the observed shape of 
the cosmic SFR history, keeping in mind however that the 
latter is still quite uncertain due to dust extinction, for ex- 
ample. If the SFR at «> 3 is as low as the UV observations 
without a dust correction suggest, and if the SFR at z < 1 is 
as high as the IR observations with the dust correction sug- 
gest, then a value for c» as small as 0.02 may be allowed. We 
also note that the c, parameter changes the size of galaxies, 
and it can hence be constrained by comparing simulation 
results with the observed scaling relations of galaxies (see 
lKobavashill2005l for the detafl). 

To investigate where the SFR is suppressed with HN 
feedback, we show the average SFRs at given total mass in 
Figure |8l At z = 6, the SFRs in massive galaxies are en- 
hanced with the metal-dependent cooling (dot-dashed and 
solid lines). After z — 3, the SFRs in low mass galaxies 
are suppressed by SN feedback (dot-dashed and dashed-dot- 
dot-dotted lines), but are unaffected in massive galaxies. If 
we include HN feedback (solid and dashed lines), the SFRs 
are smaller for all galaxies, independent of metal-dependent 
cooling. The stellar mass of the most massive galaxies are 
not ahered significantly (8.5, 6.1, 7.0, 5.8, 5.0 x lO^M© 
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Figure 8. (a) Cosmic star formation rates for SN feedback with 
(dashed line) and without (dash-dotted Hne) metal-dependent 
cooling, and for HN feedback with (solid line) and without 
(dashed l ine) it. Res t -fram e UV observati ons with dust correc- 
tion are: iLillv et aU lll996l), o pen c ircles; iMadau et al.l lll996h. 
squares; IConnoll v et al.l lll997h. filled circles: [ Steidel et al.l 



((19991) , crosse s: iBouwen s. Bro adhurst. &: lUingworthl (l2003abl'). 
filled triangles; llwata et al. 1120031). filled pentagon; iBunker et al.l 
||2004|'). filled hexa gon; iGiavalisco et al.l ||2004|'). open tria ngles; 
lOuchi et ai] ||2004| ). filled sq uares: ISchiminoyich et al.l (120051 '). dia- 
monds. Ha observations arc: Gallego et al. (^1995), filled diam ond; 
iPerez- Gonz alez ct al. (2003), Gronwall (1999), Brinchmann e t al.l 
(120041 '). iTresse fc Maddoxl (Il998l '). iTresse et al.l ( l2002t) increas- 
ing redshift with dust correction, fille d stars. Submillime - 
ter, radio and X-ray observations from: [ Hughes et al.l l|l998l) . 
open star; iBarger. Cowie. &: Richards! 1120001) . four-pointed stars; 
iNorman et al.l (|200J), asterisks, (b) Same as panel (a) but for 
the stellar density evolution. Observational data sources are: 



ICole et al.l ll200l|). fill ed circle; iBrinchrnann fc EllisI ll2q01]|). tr i- 
angles; ICohenI ( 20021 ) . open squa res; I Pickinson et al.l (12003^ . 
filled squares with t he error boxes; iRudnick et al.. (. 20061) ■ stars; 
iFontana et al.l (12009) . open circles. 



Figure 9. The average SFR normalized of the total mass for the 
models without feedback (dotted line), with SN feedback with- 
out (dash-dot-dot-dot line) and with (dot-dashed line) metal- 
dependent cooling. In addition, results are shown for HN feed- 
back without (dashed line) and with (solid line) metal-dependent 
cooling. 



respectively for the models in the label) . However, the num- 
ber of galaxies with M« ~ 10^~^^Mq decreases by a factor 
of 2 — 3, and the number of dwarf galaxies with M, < 10^ M© 
increases by a factor of 2 — 3. 

The feedback and the metal-dependent cooling affect 
the distributions of heavy elements. Figure \TU\ shows the 
redshift evolution of the metal fraction, normalized by the 
total projected metals, in all gas, in hot gas (T > lO'^ K), 
and in cold (T < 1.5 x 10^ K) gas. Because of star for- 
mation, the metal fraction increases in stars (panel a) and 
decreases in gas towards lower redshifts. At the beginning, 
metals are ejected into hot gas (panel b), which eventually 
cools and increases the metal fraction in cold gas toward 
z '^ 5 (panel c). From z ~ 5 onward, the metal fraction in 
cold gas decreases due to star formation. At z ~ 3, almost 
half of all metals are in stars. With the metal-dependent 
cooling (long-dashed and solid lines), more metals are locked 
into stars than is the case without it (dot-dashed and short- 
dashed lines). With HN feedback (short-dashed and solid 
lines), two times more metals exist in gas, more in hot 
gas and less in cold gas. This metal fraction in hot gas, 
however, is still smaller than expected from observ ations of 
the X - ray gas in cl usters (larger than two-thirds, iRenzinil 
I2OO2I ). IPave fc Oppe nhcimcr (2006) showed the increase of 
the metal fraction in hot gas toward low-redshifts using a hy- 
drodynamical model with a momentum driven wind. Such 
increase is not seen in our simulations with metal-dependent 
cooling. 

Figure [11] shows the redshift evolution of the fractions 
of metal-free gas and metal-free stars. In our model with 
HN feedback and metal-dependent cooling (solid line), half 
of the gas remains pristine until z ^ 3, and ~ 20% gas is still 
metal-free at z = (panel a) . The fraction of metal free stars 
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Figure 10. The redshift evolution of the metal fraction, nor- 
malized by the total projected metals, in (a) all the gas, (b) hot 
(T > 10^ K) gas, and (c) cold (T < 1.5 X 10* K) gas. See Fig. |8] 
for the models. 



is as large as ~ 50% at 2 '^ 6 and decreases to ~ 5% at 2 = 
(panel b). Higher resolution (short-dashed line) gives almost 
the same fraction of metal-free gas. However, the amount of 
metal-free gas strongly depends on the number of feedback 
neighbours A^fb in our model. A larger value for A'^fb gives 
a smaller fraction of metal- free gas (long-dashed line). We 
should note that we do not include an explicit mixing of 
metals between neighbouring gas particles. Gas particles are 
enriched only when they pass near dying star particles in our 
model. Although an additional mixing effect could increase 
the enriched fraction of gas, we note that the region enriched 
by a supernova is already large in our simulation model. 



4.2 Star Formation, galactic winds, and chemical 
enrichment 

In the following, we discuss in more detail the results of the 
higher resolution simulation (A'' = 96"^,A''fb = 405) with 
HN feedback and metal-dependent cooling. Figure [121 shows 
V-band images for our galaxies identified by the FOF al- 
gorithm at 2 = 0. The most massive galaxy shown (up- 
per left panel) is a merging system, with the total mass 
Mtot ~ 1.5 X 10"Mo and the stellar mass M, ~ 4 x lO^M©. 




-3 



Figure 11. The fractions of metal-free (a) gas and (b) stars. 
The solid, short-dashed, and long-dashed lines are for the mod- 
els with HN feedback and metal-dependent cooling, respectively. 
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Figure 12. V-band images of selected galaxies in our cosmolog- 
ical simulation at z = 0. The panels are 200 kpc on a side, and 
the galaxy images are ordered by their total mass. 



The others illustrate the typical morphologies of our simu- 
lated galaxies with Mtot ~ 6x 10^^ 8x lO", and 2x lO^M©, 
and M. ~ 1 X 10",4 x 10^", and 4 x IO^Mq, respectively. 
When and where do stars form? To answer this ques- 
tion, we break up the cosmic SFR history according to 
galaxy mass. In Figure 1131 we show the SFR history of sys- 
tems with final stellar masses of A/» ~ 10^^ (solid line), 
10^" (dashed line), 10^ (dot-dashed line), and IO^Mq (dot- 
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Figure 13. The cosmic SFR liistory split by different galaxy 
mass scales at 2 = 0. We show results for stellar masses of lO-"^^ 
(solid line), 10^° (dashed line), 10^ (dot-dashed Une), and 10* Mq 
(dot-dot-dot-dashed line). The dotted line shows the total. The 
galaxies are identified by EOF at z = 0, so these SFRs correspond 
to the stellar age distributions in the galaxies. 




Figure 14. The same as Fig. 1131 but the galaxies are here iden- 
tified at each redshift, so the mass quoted corresponds to the 
current one at each redshift. 



dot-dot-dashed line). The dotted line shows the total, cor- 
responding to the ~ 1000 galaxies in this simulation. The 
galaxies have been identified by FOF at 2 = for this plot, 
and thus these SFRs correspond to the age distribution of 
stars in the galaxies. For all galaxy masses, the SFRs show 
a peak around 2; ~ 3 — 4, and the majority of stars are as old 
as ~ 10 Gyr. On the other hand, in Figure 1141 we identify 
galaxies at each redshift and split up the SFRs according 
to the current stellar mass measured at the redshift, which 
are comparable to the observations of high redshift galaxies. 
This shows that most stars have formed in low-mass galaxies 
with IO^^'^Mq at high redshift 2 >3. 10^° Mq galaxies exist 
at high redshift z < 5, but 10^^ M© galaxies appear only 
after 2^2. From these two figures we conclude that most 
stars have formed in dwarf galaxies before they merge to 
massive galaxies in our simulation. As a result of the hierar- 
chical clustering of dark matter halos, such old stars belong 
to massive galaxies at low redshifts. 



If we consider luminosity-weighted ages, the largest 
galaxies need not always be the oldest if young stars keep 
on forming in them. At the present epoch, we find no re- 
lation between the luminosity-weighted ages and the stel- 
lar mass of the galaxies, and most galaxies are as old as 
~ 10 Gyr. Large galaxies with M* > 10^^ M© tend to have 
young ages as ~ 8 Gyr in our model. In these galaxies, 
feedback from active galactic nuclei (AGN) may stop re- 
cent star formation. For dwarf galaxies, the ages span a 
wide range of 1 — 10 Gyr, and there exist two popula- 
tions, young and old dwarfs that correspond to the observed 
dwarf irregulars and spheroidals, respectively. However, the 
number of young dwar fs seems to be too sm all compared 
with observations (e.g.. lKauffmann et al.ll2003r ). and the so- 
called 'down-sizing' effect of cosmic star formation is not 
clearly seen in our simulations. This is because i) star for- 
mation is not terminated in massive galaxies, presumably 
due to the lack of efficient feedback and ii) not enough 
young dwarfs are forming at 2 = 0. AGN feedback (e.g. , 
iDi Matteo. Springel. fc HernguistlboOSl . ICroton et al.ll2006h 
may solve the former problem, but not the latter. 

How are heavy elements ejected from galaxies to the 
IGM? In the simulation, we can trace the orbit of gas par- 
ticles over time. Exploiting this, we define as wind particles 
those that are not in galaxies now, but have been in galax- 
ies before. Figure [T^ shows the spatial distributions of this 
'wind gas'. The blue and red points mark the positions of 
gas particles in galaxies and of gas in the wind, respectively. 
The corresponding distribution of star particles is shown in 
the right panel. The wind particles surround the stars in 
galaxies, but do not extend very far from them. In this sim- 
ulation, ~ 10% of baryons turn into stars, ~ 10% of the gas 
stays in galaxies (~ 8% is hot), and ~ 20% is ejected as 
galactic winds (see Table [2] for the other models). The rest, 
half of the baryons, never accretes onto galaxies. 

When we follow the orbits of gas particles, we can 
also examine from which galaxies the wind gas particles are 
ejected. This allows a measurement of the ejected wind mass 
from each galaxy. It is possible that a wind gas particle is 
ejected from two different galaxies (at different times), in 
which case we count it for both galaxies. But if a wind gas 
particle should be ejected twice from the same galaxy, it 
is only counted once. In Figure 1161 we plot the wind frac- 
tion, i.e. the ratio between the total wind mass to the total 
accreted baryon mass, against the total mass. Here we plot 
only galaxies with M* > 10* Mq because of the limited reso- 
lution. Massive galaxies are dominated by baryons (6 ~ 0.1), 
but the baryon fraction is smaller for lower-mass galaxies 
(panel a), which implies that baryons are more efficiently 
ejected from them. Most of the gas in our massive galaxies 
is in the hot phase. Therefore, the difference in the baryon 
fraction does not confiict with the observed baryonic TuUy - 
Fisher relation (McGaugh 2005; Pfcnnigcr fc Revaj [JOOsI ) . 
as we explicitly show in Figure [TTl 

A clear relation is found between the wind fraction 
and the total mass (panel b). Winds are efficiently ejected 
from small galaxies, with ~ 80% of accreted baryons be- 
ing ejected from Mtot ~ 10^^ M© galaxies. A similar rela- 
tion is also found for the ejected metal fraction, i.e. the ra- 
tio between the wind metal mass to the total metal mass 
(panel c). It is interesting that the wind fraction and the 
ejected metal fraction correlate well with the stellar metal- 
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Figure 15. Spatial distribution of wind particles at 2: = 0. The blue and red points show gas particles in galaxies and in the wind phase, 
respectively. The distribution of star particles is given in the right panel, with different colours for different galaxies. 



licity (panel d). Based on this finding, we conclude that the 
origin of the mass-metallicity relation can be explained with 
galactic winds. 

We note that the mass contribution of the winds is still 
larger for massive galaxies. ~ 50% of the total wind mass 
and '~^ 40% of the total metals are ejected from massive 
galaxies with Mtot > lO^'^ Mq. The smaller contribution for 
the metal ejection from massive galaxies is due to the metal- 
dependent cooling; enriched gas tends to be locked in stellar 
population of galaxies. 

The metal enrichment timescale depends on the envi- 
ronment. Figure [TS] shows the evolution of iron and oxy- 
gen abundances in the gas phase. The average metallicity 
of the universe (solid line) reaches [0/H] ~ —2 and [Fe/H] 
~ —2.5 at z ~ 4 (1.5 Gyr), but reaches the same values at 
2 ~ 3 (2.1 Gyr) in the IGM (dotted line). The box shows 
[C/H] obtained from quasar absorption line measurements, 
which involves some uncertainty in the UV background ra- 
diation. [C/Fe] is also uncertain, and spans a wide range of 
—0.6 < [C/Fe] < 1 in the metal-poor stars in the Milky Way. 
Since core-collapse supernovae cannot produce so much car- 
bon relative to Fe ([C/Fe] ~ -0.4. iKobavashi et al.ll2006l ). 
the C-rich stars seem to be enriched by massive and/or low- 
mass stellar winds, they could also be affected by internal 
mixing or by external enrichment from a binary companion. 
The metallicity of the IGM seems not to be affected sig- 
nificantly by stellar winds because of the lower energy and 
the different timescale. Assuming that the IGM is enriched 
only by core-collapse supernovae, we could account for this 
by shifting the observed region assuming [C/Fe]= —0.5 and 
[C/0]= —1, and then our result becomes comparable to the 
observations. We note that t his assumption can explain the 
observed [Si/C] in the IGM (JAguirre et al.ll2004l '). while pair 
instability supernovae, which produce much more iron and 
less carbon, cannot resolve this confiict. A self-consistent 
model for the carbon abundance is required. 

The metallicities of the individual galaxies with Afg + 
M, > 10^ Mq are shown as points, with larger points rep- 



resenting more massive galaxies. In large galaxies, the en- 
richment takes place so quickly that [0/H] reaches ~ — 1 
at 2 ~ 7. This is consis tent with the meta llicities of LBGs 
(large errorbar at 2 = 3. |Pettini et al.ll200H ). Smaller galax- 
ies evolve with more variable timescales, resulting in a larger 
scatter of —2 < [0/H] < at present. The metallicities 
of DLA systems (err orbars, IProchaska et al.ll2003l : triangles, 
iKulkarni et al.ll2005l ) are comparable with our galaxies, sug- 
gesting an identification of DLAs with dwarf galaxies or the 
outskirts of massive galaxies. 

The metallicities measured for lo cal emission-line galax- 
ies (two parallel lines at 2 = — 1, iKobulnickv fc KewlevI 
2004) are much larger than our simulation predictions. One 
reason could be an aperture effect. The metallicities of our 
model galaxies are measured for the whole galaxy, which 
may be systematically low compared to the observations 
since strong metallicity gradients are found in our galaxies 
(Fig. [5} and the observational data focuses on star-forming 
regions which may be biased towards galaxy centres. The 
other reason is that a lot of metal-poor hot/warm gas exists 
in our large galaxies because of their deep potential. Observ- 
ing the metallicity of warm gas is also important to put con- 
straint on the chemical enrichment history of the universe. 
If we measure the metallicity of the cold gas in the inner- 
most 10 kpc region, [O/H] reaches ~ 0.2, weakly increasing 
toward 2 = 0, which is comparable to the observations. 

In Figure 1191 the metallicities of hot (dotted-open cir- 
cles) and cold (open circles) gas are shown against the stellar 
mass of the galaxies. There seems to be a relation that gives 
a larger metallicity of the cold gas in more massive galaxies, 
and [O/H] '^ and —1.5 in giant and dwarf galaxies, re- 
spectively. The filled circles show the cold gas metallicities 
within 10 kpc, which are compar able to those fou nd in SDSS 
emission- line galaxies (solid line. iTremonti et al.i ,2004'). but 
the slope is shallower. In our model, mass-metallicity rela- 
tions for central cold gas are found after 2^5, although 
there is a significant scatter. At 2 = 2, the gas metallicity is 
consistent with the observation (solid line. lErb et al.l 12006! ) 
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Figure 16. Baryon fraction (panel a), wind fraction /w = Mw/(Mg + Af, + Mw) (panel b), and ejected metal fraction /w,z 
A/z,w/(A/z,g + A/z_, + Mz^vj) (panel c) against total mass. Panel (d) shows the stellar metallicity against the ejected metal fraction. 




Figure 17. Baryonic TuU y-Fisher relation. The solid line shows 
the observation taken from iMcGauehl 1 120051 ). 



for massive galaxies, but seem to be larger than observed for 
low-mass galaxies. 

In contrast to this large scatter, tight relations are found 
between the stellar metallicity and stellar mass at any red- 
shift, as shown in Figure 1201 based on V-band luminosity- 
weighted metallicities. The mean stellar metallicity of the 
whole galaxy (open circles) reaches solar for massive galax- 
ies, and half solar for intermediate-mass galaxies. This is 
consistent with observations of nearby ellipticals when the 
effec t of metallicity gradients is taken into account (solid 
line, iKobavashi fc A rimoto 1999). The central metallicity 
measured within 10 k pc (filled circles) is two times higher 
(short-dashed line, IPahre. Diorgovski fc de Carvalhd 119981 ) 
than a global average. In our simulation, progressively more 
material has been blown away via galactic winds from 
present-day lower mass galaxies (Fig. I16p . and we find a 
mass-met allicity relation that continues for dwarf galaxies. 
However, iGallazzi et al.l ((20051) (long-dashed line) found a 
change of the slope at M* -- 10^°Mq in the SDSS, which 
does not appear in our simulations. The change of the slope 
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Figure 18. Redshift evolution of iron and oxygen abundances of 
the gas. The points show individual galaxies with the size of the 
symbols encoding the size of the galaxies. The solid, dashed, and 
dotted lines are for the mean, the largest galaxy, and the IGM, 
respectively. The obs ervational da t a are : The large errorbar at 
2 = 3 for th e LBGs oflPettirii et al.| ll200lh . errorbars for the DLA 
systems of IProchaska et al.l l|2003l ) where fO/Fe]=0.4 a dopted, 
triangles for the DLA systems of iKulkarni et al.l BOOST) where 
[Zn/Fe]=0 adopted, and the two parallel lines at z = — 1 for 
the emission-line galaxies ofjKobulnickv & Kewley 1120041 ). The 
box is for [C/H] in the IGM l|Schave et al. 2003|), shifted with our 
assumption of [C/Fe]= -0.5 and [C/0]= -1. 



Figure 19. Mean metallicities of hot (T > 10® K, dotted-open 
circles) and cold gas (T < 10'* K, open circles), plotted against the 
total stellar mass. The filled circles are for the cold gas within 10 
kpc. The solid lines show the relation for th e nebular gas observed 
with the emission lines for SDSS g alaxies iTrem onti et al.ll2004) 
at 2 = and star-forming galaxies I IErb et al.ll2006h at 2 = 2. 
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Figure 20. Relation between the V-band luminosity-weighted 
stellar metallicity and the total stellar mass. The open and 
filled circles are for the whole galaxy and central region within 
10 kpc, respectively. The solid and short-dashed lines are for 
the observed relations in nearby ellipticals with measured Mg2 
index, the average is derived accoun ting for the metallicity 
adients JK obayashi fc Arimot d 11993) and the central values 



Pahre. Diorgovski fc de Carvalhq 1998|). respectivelv . The long- 



dashed line is for SDSS galaxies IjGallazzi et al.ll200. 
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at Mt ~ 10^ Mq in our results is presumably due to resolu- 
tion effects. 



5 CONCLUSIONS 

We have simulated cosmic chemical evolution with a tree- 
SPH code including star formation, supernova feedback, and 
a detailed chemical enrichment model that does not rely 
on the instantaneous recycling approximation. We have also 
investigated the role hypernova may play for regulating star 
formation and for enriching the universe. Our main findings 
can be summarized as follows. 

(i) Supernova feedback can suppress star formation, but 
this reduction is largely compensated by an increase of sim- 
ilar size when metal- dependent cooling is included. In sim- 
ulations of isolated disk galaxies with a total halo mass of 
10^" h~^yiQ, a bipolar galactic wind is formed as a result of 
the supernovae feedback and ejects up to half of the baryons. 
The ejected metal fraction is not as large because of metal- 
dependent cooling, and amounts to a few percent for su- 
pernova feedback. If hypernovae are included as well, this 
raises to 25%. In this case, more metals are ejected with 
more energy, and the enriched gas is heated so strongly that 
rapid subsequent cooling in prevented. The wind efficiency 
does not depend strongly on numerical resolution. However, 
it depends strongly on halo mass; for massive halos of mass 
lO^'^ /i~^Mq no winds occur. 

(ii) The cosmic SFR shows a peak at 2 ~ 4, with ~ 10% of 
baryons turning into stars, roughly consistent with observa- 
tional estimates. We analyse the contributions to the cosmic 
SFR from different galaxy mass scales at present (Fig. I13p , 
and at the observed redshift (Fig. I14[) . The majority of the 
stars in present-day massive galaxies have formed in smaller 
systems at high redshifts, and the stellar population of these 
galaxies is as old as 10 Gyr, consistent with observations of 
elliptical galaxies. No relation is found between luminosity- 
weighted age and stellar mass. The ages of low-mass galaxies 
span a wide range of 1 — 10 Gyr. 

(iii) Tracing the orbits of gas particles, we found that 
~ 20% of all baryons are ejected at least once from galaxies 
into the IGM. Galactic winds blow particularly efficiently 
from low-mass galaxies. The wind fraction is larger for less 
massive galaxies, and ~ 80% of the accreted baryons are 
ejected from galaxies with mass below Mtot ~ 10^^ M© 
(Fig. I16|l . The wind fraction correlates well with the stellar 
metallicity, and the origin of the mass-metallicity relation is 
from galactic winds. 

(iv) The metal-enrichment history depends on environ- 
ment. In large galaxies, enrichment takes place so quickly 
that [0/H] reaches ~ — 1 at z ~ 7, which is consistent with 
the sub-solar metallicities of the Lyman break galaxies. The 
low metallicities of DLA systems are also consistent with 
our galaxies, provided these systems are dwarf galaxies or 
the outskirts of massive galaxies. The low [C/H] of the IGM 
can be explained if the IGM is enriched only by SNe II and 
HNe. The average metallicity of the universe reaches [O/H] 

-- -2 and [Fe/H] 2.5 at 2 ~ 4 (1.5 Gyr), but reaches 

the same values at 2 ~ 3 (2.1 Gyr) in the IGM. 

(v) In galaxies, metallicity of the cold gas increases with 
galaxy mass, which is comparable to observations with a 
large scatter. The central cold gas shows a relation between 



galaxy mass and metallicity with shallower slope than ob- 
served in emission-line galaxies. For the stellar population, 
the observed mass-metallicity relation is well reproduced, 
and originates in mass-dependent galactic winds. These re- 
lations are present since z ~ 5. 

At z = 0, however, star formation has not terminated 
in our massive galaxies, and not many dwarf galaxies are 
still forming. The latter problem may be due to lack of reso- 
lution, while the former may reflect non-supernova feedback 
such as accretion onto supermassive black holes. It will be 
interesting to include active galactic nuclei in future studies 
of cosmic chemical evolution, as well as to further increase 
the dynamic range of these calculations. This will then allow 
methods such as the ones explored here to fully exploit the 
observational abundance measurements to better constrain 
galaxy formation physics. 
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